Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Spectral analysis of the sdO K648, the exciting star of the 
planetary nebula Ps 1 in the globular cluster M 1 5 (NGC 7078)^ 

T. Rauch\ U. Heber**^, and K. Werner^ 

^ Institut fiir Astronomie und Astrophysik, Universitat Tubingen, D-72076 Tubingen, Germany 
^ Dr. Remeis-Sternwarte, Universitat Erlangen-Nurnberg, D-96049 Bamberg, Germany 

Received 17 July 2001 / Accepted 26 October 2001 

Abstract. We present a spectral analysis of the sdO central star K 648 based on medium-resolution optical and high-resolution 
UV spectra. The photospheric parameters are determined by means of state-of-the-art NLTE model atmosphere techniques. We 
found Tcft — 39±2kK and log 17 = 3. 9 ±0.2. The helium (nHe/nn ~ 0.08) and oxygen (no/nu = 0.001) abundances are about 
solar while carbon is enriched by a factor of 2.5 (nc/jiH = 0.001). Nitrogen (nN/wn ~ 1 ■ 10"^, [N/H] = -2.0) appears at a 
sub-solar value. However, these metal abundances are much higher than the cluster's metallicity (M 15: [Fe/H] = -2.25). The 
surface composition appears to be a mixture of the original hydrogen-rich material and products of helium burning (3a process) 
which have been mixed up to the surface. The abundances of He, C, and N are consistent with the nebular abundance, while O is 
considerably more abundant in the photosphere than in the nebula. From a comparison of its position in the log Tcff-log g plane 
with evolutionary calculations a mass of O.STIq qi ^0 ™'l ^ luminosity of 3810±1200Lq are deduced. Our spectroscopic 
distance d = ll.lj^j g kpc is in agreement with the distance of M 15 as determined by Alves etal. (2000). From the GHRS 
spectra we measure a radial velocity of Vmxi = —130 km/sec. 

Key words. Galaxy: globular clusters: individual: M 15 - ISM: planetary nebulae: individual: Ps 1 - Stars: abundances - 
Stars: AGB and post-AGB - Stars: evolution - Stars: individual: K 648 



1. Introduction 

K648 (my = 14.73, Alves etal. 2000) was registered by 
. Kustner (1921) close to the center (Aa = +14", AS = 
' +26") of the galactic globular cluster M 15. A spectrum 
taken by Pease (1928) shows a continuous O-type spectrum 
and characteristic emission lines of a planetary nebula (Ps 1, 
PNG065.0-27.3). Pease derived a radial velocity of t;rad ~ 
— 180±50km/sec from his spectrum and concluded that K 648 
is probably the first PN discovered in a globular cluster 

Joy (1949) improved the radial-velocity measurement 
(viad = — 115km/sec) and concluded that, without doubt, Ps 1 
is a member of M 15 (Wrad — — 107km/sec, HaiTis 1996) be- 
cause of the close agreement of their radial velocities. 

The existence of a PN in M 15 is unexpected: At a turn-off 
age of 12 billion years, the most massive main-sequence stars 
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in M 15 should have Afinitiai~0.8 Mq. For such low-mass stars 
it is almost impossible to ascend to the AGB and eject a PN. An 
analysis of the nebula properties by Adams et al. (1984) yields 
a solar He/H ratio and a slightly higher than solar C/H ratio 
while the total N, O, and Ne abundance is less than solar by 
a factor of 18 (Fig.|l3|). The high carbon abundance in the PN 
was interpreted by Adams etal. as a product of helium burn- 
ing (3q; process) and a subsequent third dredge-up which has 
brought the material to the stellar surface. The nebula ejection 
has taken place then after this event. 

Alves etal. (2000) suggest that the progenitor of K648 
experienced mass augmentation in a close binary merger and 
thus, the remnant has a higher mass than remnants of single 
stars in M 15. 

A preliminary analysis of the optical and GHRS spec- 
tra (Tab. ^ by means of line-blanketed NLTE model atmo- 
sphere techniques by Heber et al. (1993) revealed Tcff = 37kK, 
log g = 4.0, riHe/^-H ~ 0.5 (by number) and a three times solar 
carbon abundance. Since this is much higher than the metal- 
licity of M 15 an explanation for the evolutionary history of 
K648 and its PN is difficult. Bianchi etal. (1995, see this pa- 
per for a more detailed introduction) presented a first anal- 
ysis of HST data and aiTived at Tcff = 35kK. Another spec- 
tral analysis based on high-resolution Keck spectra was pre- 
sented by McCarthy et al. (1997). They arrived at Tcff = 43kK, 
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Fig. 1. Optical spectrum of K648 obtained with the TWIN spectrograph (top). Positions of H I, He I-II, and C III-IV lines are 
indicated. Note the contamination by nebula emission, e.g. around H 7 and He i A 447lA. A synthetic spectrum (parameters from 
Tab. ||) is shown at the bottom 



log 5 = 3.9, nHo/"-H = 0.08 and a solar carbon abundance. 
However, in contrast to Heber etal. (1993), NLTE line blan- 
keting had not been taken in account in these models. 

Recently a LTE abundance analysis was presented by 
Bianchi etal. (2001): They used Tog and \ogg from Heber 
et al.(1993) and arrived even at a two times higher helium abun- 
dance of nHo/'T-H = 0.6, carbon is four times solar, oxygen un- 
derabundant, and silicon is about solar 

In order to investigate the enigma of K648, we have per- 
formed a new NLTE spectral analysis based on state-of-the-art 
metal- line blanketed model atmospheres. 



and the interstellar reddening (Sect. jzO ). The NLTE spectral 
analysis is then described in Sect.|3| 



2.1 .1 . Interstellar sulfur lines 

Two interstellar lines, S 11 AA 1250.59, 1253. 8IA, are promi- 
nent in the GHRS spectra. This allows to measure the wave- 
length shift of the ISM: It is blueshifted by about 0.145 A (this 
is equivalent to Uioi = — 35km/sec). This values allows e.g. 
to identify the interstellar contribution to the photospheric N v 
lines (Sect.O). 



2. Observations 

2. 1. Ultraviolet spectra 

The following spectral analysis is mainly based on HST "pre- 
view" spectra (Tab. |l]) which were retrieved from the archive 
operated at ST-ECF. 

The high-resolution GHRS spectra allow to measure 
the radial velocity of K648 precisely: In Z1BN0106T and 
Z1BN0107T (Tab. the photospheric lines are blueshifted by 
0.55 and 0.57A (we used the C iii C iv, N v, and O iv lines in 
the spectra to measure the shift) and hence Vmd = —128 and 
— 133 km/sec, respectively. This is close to the values given 
by Joy (1949, -115 and -129 km/sec) and Schneider etal 
(1983, —128 km/sec). However, the velocity of M 15 is some- 
what smaller (Vrad = —107 km/sec, Harris 1996). 

The analysis of UV spectra requires a careful determination 
of the stellar continuum. Thus we start already here w ith an 
investigation of the interstellar H i column density (Sect. 2.1.2) 



2.1.2. Interstellar neutral hydrogen and reddening 

A Hi column density of nni = ^(^1) ' lO^^cm^^ is de- 
termined from the Lya profile (Fig.g). The interstellar red- 
dening is then estimated by means of the relation nni = 
(3.8±0.9)-1021£;b_v (Groenewegen & Lamers 1989) yielding 
a color excess of Eb-v — 0. 13lg Q4. We use this result to com- 
pare a synthetic spectrum (from a H-nHe model) with the FOS 
spectrum of K648 (Fig.|]). In agreement with the above val- 
ues, we achieve a good fit at Eb-v — 0.10 in the wavelength 
range 2100-3200 A. This is also the amount of the foreground 
reddening towards M 15 (Harris 1996). 

At shorter wavelengths the model predicts a higher flux 
level than observed. This discrepancy is not caused by metal- 
line blanketing. Even when we increase the abundances to so- 
lar, no match can be achieved (Fig. Q). Instead there a signifi- 
cantly different reddening law might be valid. 

Since lUE spectra (e.g. SWP 17069) show the same flux 
level like the FOS spectra in the wavelength range around 
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Fig. 2. Comparison of the theoretical Ly a profile with the FOS 
spectrum. The best fit is achieved at nni = 5 ■ lO^^/cm^. The 
theoretical flux is reddened according to rt H i (th e continuum is 
scaled to fit the line wings of Ly a. Sect. 2.1.2 ). The parame- 
ters of the line-blanketed model atmosphere are summarized in 
Tab. I 

1400A (not shown), problems in the flux calibration appear to 
be unlikely. 

In the FOS spectra, we will use the He ii Fowler series for 
our analysis (Fig.^ which is located in the part of the spectrum 
that is well matched by the model. Thus we adopt E'b-v — 
0.10 for our analysis. In the case of detailed line profile fits 
to the FOS and GHRS spectra, the model continuum has to 
be individually scaled in order to fit the observed continuum 
around the analyzed line. 

2.2. Optical spectra 

The optical spectra (resolution = 1 A) of K 648 were obtained 
at the 3.5m telescope at the Calar Alto (DSAZ, Spain) on June 
12, 1989 with an exposure time of 2100 sec. Since these spec- 
tra are, in contrast to the UV spectra which are used here, not 
absolutely calibrated we use the rectified spectrum (e.g. Fig.[l|) 
instead. 

Since K 648 lies close to the center of M 15, the background 
determination was difficult (Heber etal. 1993) due to the low 
spatial resolution of the TWIN spectrograph. 

Some of the "strategic" spectral lines are contaminated 
by nebular emission (Fig.|l|) and thus unreliable: E.g. the 
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Fig. 3. Determination of the interstellar reddening. A syn- 
thetic spectrum, calculated from a H-nHe model (Tcff = 39kK, 
log (7 = 3.9, riHc/^-H = 0.08), is shown with a reddening ac- 
cording to -Eb-v = 0.07, 0.10, and 0.13. At longer wavelengths 
(A~2000A) the observation is well fitted with Eb-v = 0.10 
while at shorter wavelengths the reddening is underestimated. 
In the optical wavelength range (bottom), the measurements by 
Alves et al.(2000) are indicated with their error bars. The higher 
i?B-v seems to fit slightly better in the optical, however, then 
the fit to the measured UV flux is worse 



He I / He ii ionization equilibrium which is usually a very sensi- 
tive indicator for the effective temperature Tcff cannot be eval- 
uated because all He i lines are filled in by nebular emission. 
However, He ii lines as well as some carbon lines can be mea- 
sured and will be used in our spectral analysis. 



3. NLTE analysis 

Bianchi et al. (1995) found evidence for a stellar wind of K648 
and derived Vinf = 1630km/sec and log Af /(Moyr^^) = 
—8.2 from the strong C iv AA1548,1550A resonance doublet. 
Other spectral lines appear unaffected by the wind. Hence, the 
plane-parallel, hydrostatic NLTE model atmospheres (Werner 
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Fig. 5. UV spectrum of K648 obtained with the FOS (top). Same Uke Fig.|l| Note that the Hen Hnes, calculated with 
nuc/nn = 0.08 are well matched 



1986, Rauch 1997, Rauch 2000, and references therein) are ap- 
propriate for this analysis. 

In order to model the background opacity, metal-line blan- 
keting of the iron-group elements (Dreizler & Werner 1993, 
Haas et al. 1996) is considered with a generic model atom. This 
is constructed from all elements Ca - Ni (ionization stages III - 
VII with all lines given by Kurucz 1996). Reduced abundances 
([Z/H] = -2.25, following Harris 1996) are assumed. However, 
the impact of the Ca-Ni lines on the spectrum is small (Fig. 

In order to check the validity of the preliminary pa- 
rameters from Heber etal. 1993 (Toff = 37kK, log g = 4.0, 
JT-Hc/^-H = 0.5) we have performed some test calculations. 
From the comparison of synthetic spectra calculated from the 
new H-nHe models to the observation we find that g is in good 
agreement, while T^ff is higher and riHc/"-H is about solar. The 
differences can be explained by the fact that our NLTE model 
atmospheres have been significantly improved since 1993 (e.g. 
Rauch 1997, Rauch 2000): The metal-line blanketing of all el- 
ements up to the iron group is considered in detail and thus 
the temperature structure is correctly modelled and subsequent 
line-formation calculations in order to determine metal abun- 
dances are much more reliable. E.g. the high He abundance 
found by Heber et al. (1993) and Bianchi et al. (2000) is due to 



artificial effects if model atoms are too small and/or metal-line 
blanketing is not accounted for in the model atmosphere calcu- 
lation. 

Due to the impact of CNO on the atmospheric struc- 
ture and background opacity we decided to calculate a grid 
of HH-HeH-CH-NH-OH-(Ca-Ni) models in order to use e.g. the 
C iii/C IV ionization equilibrium to determine T^e reliably. The 
analysis is described in the following. 

3. 1. Surface gravity and helium abundance 

The first parameter to be fixed is g. For this purpose we 
use two strongest uncontaminated lines in the optical spec- 
trum. He II AA 4200, 454lA (Fig. |l]). We have calculated H-nHe 
models with Teg = 37 - 43kK, logg = 3.5-4.5 at different 
nuc/nn ratios. Although Teg cannot reliably be judged from 
this approach, the riHo/^-H is about solar. He ii A 4686A may be 
contaminated by nebular emission is not considered. Models 
with log g = 3.7 — 4.2 and nHo/"-H = 0.06 — 0.10 (by num- 
ber) fit the observation well (Fig.^). These values are con- 
sistent with our more detailed models which account for 
metal-line blanketing (see below). We adopt logsr = 3.9 and 
"-Ho/"-H = 0.08. An error of 0.3 dex can be estimated by the 
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Fig. 4. Impact of iron-group opacities on the UV flux (model 
parameters: Tab. Even at solar iron-group abundances (thick 
line), the model flux is much higher than observed. At the clus- 
ter metallicity [Fe/H] = -2.25, the iron-group opacities can be 
almost neglected (fine line). The synthetic spectrum is scaled 
to fit the observed flux at 3 000 A 

variation of Tcff and log g within their error range and the qual- 
ity of the available spectrum. 

3.2. Effective temperature and carbon abundance 

In the optical spectrum of K648, the He i lines are contami- 
nated by nebular emission, and thus, a commonly used tool to 
determine Tcff, the He i/He ii ionization equilibrium, cannot 
be used. The GHRS UV observations (Tab.[l|) had been used 
with the aim to evaluate the C iii / C iv ionization equilibrium 
(Heber etal. 1993) which is also a very sensitive indicator for 
Teff(cf. Rauch 1993). 

Since there are a lot of carbon lines identified in the op- 
tical spectrum (Fig.|l]) we will use them to determine the car- 
bon abundance first and then model the C iii / C iv ionization 
equilibrium precisely. New HH-HeH-CH-NH-OH-(Ca-Ni) models 
(Toff = 39kK, log 3 = 3.9, nHc/^H = 0.08, [(Ca-Ni)/H] = -2.25 
have been calculated and the CNO abundances were adjusted to 
fit the observation (Fig. We adopt ric/nn = 0.001 ±0.3 dex. 

We calculated synthetic spectra from a small grid of models 
with Teff = 37 - 43kK and nc /^h = 0.001. It turned out that 
the C in lines are strongly dependent on variation of Tcff while 
the C IV lines appear almost unchanged (Fig.^. From the fit 
of the strongest C iii line in the FOS spectrum of K648, the 
C III A 1175A triplet, we determine roff = 39±2kK (Fig.||). 

This result is verified by the C iii and C iv lines in the high- 
resolution GHRS spectra (Fig.^. 

3.3. Photospheric nitrogen and oxygen abundances 

In the last step, we determine the N and O abundances. New 
HH-HeH-CH-N-nO + iron-group (Sect.||) models are calculated. 
In the GHRS spectra, the most prominent N and O lines are the 
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Fig. 6. Theoretical line profiles (Toff — 39kK, convolved with 
a Gaussian of lA FWMH in order to match the instrumental 
resolution) of He ii A 4200A (top) and He ii A 454lA (bottom) 
compared with the observation. At 3.7~log(7 or log5~4.2 no 
good fit can be achieved 

Nv A1238,1242A resonance doublet (apparently unaffected 
by the stellar wind (Sect.|) and the O iv AA 1338, 1343A lines. 
We use these in order to determine the abundances. We achieve 
nN/nn = 1 • 10-*^ (Fig.|l§) and no/nH = l-10-3 (Fig.[l]). 
Due to the S/N ratio of the spectra, an error of ±0.5 dex has 
to be assumed. Note that the line cores are much too strong at 
higher abundances. 

4. Mass, luminosity, post-AG B age, and 
spectroscopic distance 

In Fig. |l^ we compare the position of K 648 with standard evo- 
lutionary tracks of hydrogen-burning post-AGB stars. From 
the evolutionary calculations of Schonberner (1983) we in- 
terpolate a stellar mass of 0.57t[5:!Ji Mq, a luminosity of 
3 810±1200 Lq, and a post-AGB age of 6 800l^^j^^ a. 

From our final model (Tab.|), we can determine the spec- 
troscopic distance of K 648 using the flux calibration of Heber 
etal (1984): 

/v = 3.58 • 10"^ • 10"°-^™^oergcm~2 sec"! k'^ (1) 
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with mvo = my - 2.175c, c = 1.47£;b-v = 0.147 (Fig.|), 
my — 14.73, and AI ~ O.57M0, the distance is derived from 
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changed in this Teff range 
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Fig. 10. Synthetic spectra with rcff = 39kK, log .g = 3.9, 
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Alves etal.2000) it is 0.166to^i^| x 0.145tg:o^8 PC- If we as 



d= 7.11- IQ^\Jh^ ■ M ■ 100-4™vo-log3 pc. 



(2) 



044 ^ "-^^^-0.038 J 

sume an average expansion velocity of Ve^^-p = 20 km/sec, then 
the dynamical age is tdyn — 3950^^050 ^- ^^^^ is within the er- 
ror range in agreement with its post-AGB age (6 800^2 loo ^' 
Sect.H). 



With the Eddington flux at Acff = 5454A our final model atmo- 
sphere Hi, = 5.31 • lO^^'erg cm~^ sec^^ Hz^^ we derive a dis- 
tance of d = 11.1^*^21) kpc- This distance is in agreement with 
the distance of the globular cluster M 15 (rf = 12.3±0.6kpc, 
Alves etak 2000). 

5. The surrounding nebula Ps1 

For the ambient nebula, we can calculate the following: The lin- 
ear dimension of the inner shell (0'.'8 x 0'.'6, Alves et al. 2000) is 
OMitoml X 0m2toml pc, and of the outer shell (3'.'1 x 2'.'7, 



6. Results and discussion 

In our analysis it became clear that the determination of T^s 
from the two available GHRS spectra alone (C iii / C iv ioniza- 
tion equilibrium) is not precise enough (too few lines to evalu- 
ate) - together with the available lUE, FOS, and TWIN spectra, 
however, we determined Teff = 39±2kK and log g = 3. 9 ±0.2 
which is within the error ranges in good agreement with the 
preliminary values given by Heber et al. (1993, 37 kK / 4.0). In 
contrast to Heber etal., however, we arrive at an almost solar 
("-Hc/'^H = 0.08) He abundance. The high carbon abundance 
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Fig. 12. Position of K648 in the logTcff-logg plane com- 
pared to theoretical evolutionary tracks of post-AGB hydrogen- 
burning stars (dashed, masses given in M0, Schonberner 1983) 
and post-EHB evolutionary tracks for [Fe/H] = -1.48 (Dorman 
etal. 1993). The positions of other UV-bright stars in globular 
clusters (cf. Moehler etal. 1998, Fig. 3) are marked. The cross 
in the lower right part indicates the typical errors 

(2.5 times solar) is verified. The photospheric parameters of 
K 648 are summarized in Tab. |[ 

The solar riHo/^iH surface abundance indicates that K648 
might be a H-burning post-AGB star. Compared to the metallic- 
ity of M 15, the C and O abundances appear enriched (Fig. pj] ), 
likely due to a He-shell flash. 

In order to further improve the analysis and the abun- 
dance determinations, high-resolution and high-S/N spectra 



Table 2. Parameters of K 648. The errors are estimated by the 
variation of Toff, logg and the abundance ratios within their 
eiTor limits in order to achieve a fit of the synthetic spec- 
trum to the observation. The solar abundance ratios are given 
by Holweger (1979) and Sturenburg & Holweger (1990). The 
cluster metallicity is times solar (Harris 1996) 





Teff/kK 


39.0 


±2.0 
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3.9 
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8 • 10"^ 


±0.3 dex 
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Fig. 13. Photospheric abundances of K 648 ( 


*: this work. A: 



Bianchi etal. 2001) and its PN (o: Adams etal. 1984) relative 
to the sun. The metalicities of the sun and of M 15 are indicated 



with high spatial resolution (minimize the contamination by 
the nebula) are highly desirable: Such spectra which cover 
the complete UV range (many lines of different elements and 
different ionization stages improvement of Tog, cf. Rauch 
1993) will separate the possible interstellar and photospheric 
lines, e.g. in N v AA 1238, 1242A. Moreover, we can identify 
Si IV AA 1393, 1402A and S v A 1502A in the available FOS 
spectra at a low S/N ratio which hampers a reliable analysis. 
Better spectra allow a precise abundance determination of these 
elements. The optical wavelength range (down to the B aimer 
edge — > improvement of g) gives additional constraints. 
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